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DEVELOPMENT  AND  AERODYNAMICS  Op  PROBOSCIOID  ABLATING  NOSE  TIPS 


ABSTRACT 


\y  Rationalization  of  slender  concave  shapes  on  ablating  nose  tips  is 
obtained  In  terms  of  composite  laminar/turbulent  self-preserving 
solutions  to  the  nose  recession  problem.  The  occurrence  of  sub- 
critical  and  critical  (violently  oscillatory)  flow  conditions  is 
discussed  and  associated  with  specific  geometrical  characteristics 
of  the  ablating  nose  shape.  The  presence  of  distinct  oscillatory 
nodes  Is  noted  for  both  the  subcritical  and  critical  cases,  and  the 
distinct  underlying  inviscid  instability  mechanisms  are  pointed  out. 
On  that  basis  numerical  flow  analyses  by  extension  of  available 
unsteady  Inviscid  flow  models  is  suggested,  and  the  nature  of  the 
required  .extensions  Is  outlined. 
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I.  INTRODUCTION 


Recent  laboratory  and  flight  experiments  have  exhibited  large  departures 
between  predicted  (Refers.  1 and  2)  and  observed  (Ref.  3)  equilibrium 
shapes  for  ablating  nose  tips,  such  as  are  used  on  high  performance 
ballistic  reentry  vehicles.  Since  reliable  estimates  of  nose  tip  con- 
figuration and  survival  are  essential  to  the  determination  of  vehicle 
aerodynamics,  the  divergence  between  theory  and  experiments  must  be 
resolved  and  improved  predictive  methods  of  ablative  shape  change  need 
be  developed  with  attention  to  the  novel  aspects  of  flow  field  behavior 
(e.g.,  violent  oscillations)  revealed  in  the  experiments. 

A two  step  approach  to  the  problem  has  been  followed  in  the  investigations 
reported  here.  A qualitative  theoretical  rationalization  of  the  observed 
shapes  has  been  sought  first.  Upon  resolution  of  this  matter, the  dynamics 
of  the  associated  flow  fields  have  been  examined  and  the  formulation  of 
relevant  nunerical  models  has  been  extracted.  The  first  phase  has  largely 
depended  upon  theoretical  considerations  based  on  approximate  analytical 
models.  By  contrast,  the  second  phase  has  largely  relied  on  data  analyses. 

The  lines  of  Inquiry  followed  in  the  two  phases  of  the  research  and  the 
related  results  are  presented  sequentially  in  the  main  body  of  the  report. 
First,  rationalization  of  the  slender  concave  configurations  observed  in 
the  experiments  is  obtained  on  the  basis  of  self-preserving  (equilibrium) 
solutions  to  the  nose  recession  problem  under  conditions  where  transition 
from  laminar  to  turbulent  flow  regime  occurs  on  the  nose  (Section  II). 

Upon  the  determination  of  equilibrium  solutions  which  exhibit  the  desired 
geometrical  characteristics,  wind  tunnel  data  for  typical  configurations  of 
interest  are  examined  with  specific  attention  to  the  presence  of  oscillations 
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and  the  nature  of  the  underlying  instability  mechanisms  (Section  III).  In 
this  process  the  association  of  subcritical  and  critical  (violently  oscillatory) 
flow  situations  with  specific  geometrical  characteristics  of  the  ablating 
nose  shape  is  established.  The  development  of  either  family  of  shapes  as 
a function  of  transition  Reynolds  number  vis  a vis  flight  Reynolds  number  is 
discussed.  The  essentially  inviscid  nature  of  the  instability  mechanisms 
underlying  the  oscillations  is  pointed  out.  Based  upon  this  last  determination 
requirements  for  the  numerical  analysis  of  the  considered  flow  fields  are  set 
forth;  extensions  of  available  unsteady  inviscid  flow  techniques  meeting 
those  requirements  are  outlined  (Section  IV). 


-2- 


2 


II.  SELF  PRESERVING  ABLATING  NOSE  TIPS 


Whereas  previous  theoretical  studies  of  equilibrium  (self-preserving) 
nose  shapes, In  either  strictly  laminar  (Ref.  1)  or  turbulent  (Ref.  2) 
flow  .Identified  only  blunt  configurations,  the  experiments  of  Reference  3 
have  shown  that  relatively  slender  concave  shapes  also  develop  within  particu- 
lar ranges  of  Reynolds. number.  In  fact,  two  classes  of  such  shapes  arise 
(Figure  1):  the  first, sometimes  designated  as  sub-critical  turbulent,  is 

characterized  by  relative  overall  slenderness,  and  steady,  either  attached 
or  separated,  boundary  layer  around  the  concave  portion  of  the  nose;  the  second, 
sometimes  designated  as  critical  turbulent,  is  characterized  by  relative  over- 
all bluntness,  separated  boundary  layer  around  the  concave  portion  of  the 
nose,  unsteady  shock  layer  and  random  nose  tip  asymnetries. 

An  examination  of  the  data  readily  reveals  a distinction  between  the  viscous 
flow  regimes  encountered  in  the  experiments  (Ref.  3)  and  assumed  in  the 
analyses  (Refs.  1 and  2).  In  the  experiments,  transition  from  laminar  to 
turbulent  regime  probably  occurs  near  the  beginning  of  the  concave  portion 
of  the  nose;  in  the  analyses  a single  viscous  flow  regime  is  considered.  Thus, 
a resolution  of  the  problem  is  suggested  in  terms  of  matched,  spatially 
sequenced, nose  contour  portions,  respectively  described  by  laminar  and  tur- 
bulent self-preserving  solutions.  The  possibility  is  analyzed  below.  For 
that  purpose  the  existence  of  self  preserving  laminar  and  turbulent  solutions 
consistent  with  concave  nose  shapes  Is  examined  first.  The  matching  of  a 
laminar  and  a turbulent  solution, to  yield  configurations  of  the  type  shown 
In  Figure  1,1s  then  considered.  x 


With  a view  to  obtaining  closed  form  results  the  analysis  is  carried 
out  under  several  simplifying  assumptions, namely;  . l)perfect  gas 
flow  with  hypersonic  free  stream  velocity;  2)  pressure  distribution  on  the 
body  described  by  Newtonian  theory;  3)  boundary  layer  flow  attached  to 
the  body  and  amenable  to  representation  by  local  similarity  solutions; 

4)  uniform  wall  enthalpy  and  effective  heat  of  ablation  per  unit  volume 
of  the  surface  material.  These  approximations,  almost  classical  in  studies 
of  viscous  hypersonic  flows,  and  specifically  used  in  References  1 and  2, 
generally  yield  results  of  sufficient  accuracy  for  purposes  of  problem  defi- 
nition, as  addressed  here.  With  the  notation  of  Figure  1 for  coordinate  sys- 
tem and  body  geometry,  the  local  laminar  (qt)  and  turbulent  (q^,  flat  plate 
reference  enthalpy)  heat  transfer  rates  are  then  given  by  (Refs.  4 and  2). 


• 0.33  Pr-2/3  (Hfi  - Hw)(pe/p0)(Ue/U-)rJ 


“ ^ 

/.  (pe 


P„  u U 

* O a 


1 1/2 


Wu-)r2jds 


a) 


qt  - 0.03  Pr-^3  (Hg  - Hj  (Pf  Up  s/ur)°-e(Mr/s)  (2) 

where  H,  p,  U,  p,  v , Pr  respectively  denote  stagnation  enthalpy,  pressure, 
velocity,  density,  viscosity, Prandtl  number,  subscripts  t,  t respectively 
Identify  laminar  and  turbulent  flow  situations,  subscripts  e,  w,  r,  0 respec- 
tively characterize  flow  properties  at  the  edge  of  the  boundary  layer,  at  the 
well,  at  an  intermediate  reference  condition,  and  at  the  stagnation  point. 
Pressure  and  velocity  distributions  at  the  edge  of  the  boundary  layer 
•re  described  by 

(pe/P0  ) * cos2® 

(l^/Uj  * [1  - (cose)  ji/2 


for  blunt  configurations,  and  by 


(3a) 

(3b) 


i 


(4a) 

(4b) 


(Pe/P0>  = cos20 
(Ue/Uw)  = 1 + O{c°s20) 
for  slender  configurations. 

Self-preserving  shapes  are  determined  by  solutions  of  the  equation 

(q/cose)  = i Heff  = COnstant  <*) 

with  q described  by  either  equation  (1)  or  (2)  in  conjunction  with  either 
(3a, b)  or  (4a,b). 

Substitution  of  (1)  into  (5)  yields  the  following  equation  for  the  self  pre- 
serving laminar  nose  shape. 

f2(e)  = / r f(e;  r2J  dr 


where  C is  the  constant 


Ct  « 0.33  Pr-2/3 (Hft  - Hw)p0u0U_  (cose  /qt) 


and 


f(e)  * (Pe/P0)(ue/u.)  (cos  0)‘1 


(6) 

(7a) 

(7b) 


a function  rendered  explicit  by  (3a, b)  and  (4a, b). 

The  solution  is 

r * 2(j  + 1)  Ctf(e) 

complemented  by  the  geometrical  relation 


x * / tane  dr 
o 


Upon  substitution  of  (3a, b)  and  (7b)  into  (8a, b)  the  blunt  self  preserving 
laminar  shape  is  obtained,  viz. 


(8a) 


(8b) 
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■ Rn  ^ rj^tsine  [1  - (cos©)2(y 

• / (cose)'l[l-(cos©)2^r'*^Y],*ds} 
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where 


**  ■ *<i  * » (^,W 


1*  the  (arbitrary)  radius  of  curvature  of  the  nose  at  the  stagnation  point. 
Inversely  proportional  to  the  square  of  the  recession  rate  (qt/cos  e Hftff). 
The  magnitude  of  RN  and  of  the  recession  rate  become  uniquely  defined 
only  upon  stipulation  of  the  value  e-»  e$  at  the  nose  shoulder.  In  that 
connection  It  Is  noted  that  the  blunt  self-preserving  shape  (9a, b)  exhibits 


a maximum  In  r for 


. / v \T/a(v-i) 

• - ®i  * eos  (n) 

e.g. , « 49,9*  for  V ■ 7/5.  Mlth  attention  restricted  to  values  i et  , 

the  magnitude  of  the  radius  of  curvature  Rj,  (the  recession  rate  qt/ctfs*r«eff) 
Is  minimised  (maximised)  by  the  choice  e$  » et«  A"  alternate  choice  (Refs.  1 
•nd  2)  is  to  set  6S  « 6*.  e*  being  the  body  inclination  at  the  sonic  point 
defined  by  the  assumed  Newtonian  pressure  distribution  (3e),  vis. 

» cos*lt2/(f  ♦ l)]T/l(t  * X)  (12) 

which  satisfies  the  condition  e*  « ©t  for  1 J y 'f  2.  For  y * 7/5  the 
tMO  assumptions,  either  ©s  ■ ©t  or  «s  » e*,  lead  essentially  to  the 
self -preserving  s.tape  as  shown  in  Figure  2. 

The  daest  for  slender  self-preserving  laminar  shapes,  defined  by  equations 
(4a,  b)  (7b)  and  (8a, b),  proves  to  be  fruitless  because  logarithmically 


• ■ 9. 


divergent  values  of  x are  obtained  as  r o ande  « cos-lr  -►  (»/ 2). 

Substitution  of  (2)  Into  (5)  readily  yields  an  algebraic  relation  between 
S and  e for  a self  preserving  turbulent  shape,  viz. 

a • Ct(wr/v0)C(T0/Tr)(Pe/p0)(U#/lO]“  (cose)-* 

vw»re  Ct  Is  the  constant 

Ct  • [0.03  Pr  -2/J(He  - Hw)u0(eoso/qt)]5  [p0  U./v0]fc  (13) 

With  the  further  reasonable  assumption  Tf«  TQ  the  turbulent  self-preserv- 
ing shape  Is  described  by 


s • Ceese)*  [1  - (cosa) 
r • tl  • (eosej  ■***  - 

(«K®r  j}  - g-r4  t<»«]  icoS«i'<’r  * 


( 


and  the  slender  self-preserving  turbulent  shape  (e  -*■  v/2  as  s -*■  o) 

\3  (16a) 


s 

r 


Ct  (cos  0)3 

(3/4)  Ct(cose  )*  (16b) 

x * (3/16)Ct  { sine[4  (cose)3  - 2 cose  ] + ir  - 29}  (16c) 

are  obtained,  respectively.  The  results  (15,a,b,c)  are  displayed  in..Figure  3 
for  Y = 7/5,  while  the  results  (16a, b,c)  are  shown  in  Figure  4 for  a general  Y. 
The  blunt  solution  is  only  applicable  in  the  range  0 £ e £ et  because 
(ds/de)  becomes  null  at 


cos 


-i  [ Y ± 2 

[ 7y  - 4 J 


y/2(y  -l) 


e.g. 


et  * 55.36°  for  Y = 7/5  . By  contrast,  the  slender  solution  is 


applicable  to  the  entire  range  U/2)  £ e £ 0.  As  in  the  laminar  case, 
the  magnitude  of  Ct  (which  has  the  dimensions  of  a length)and  of  the  reces- 
sion rate  remain  arbitrary  until  the  value  e = es  at  the  nose  shoulder  is 
specified.  For  the  blunt  shape  the  magnitude  of  (the  recession  rate 
q^/cose  Heff)  is  minimized  (maximized)  by  the  choice  ©s  = e^;  for  the  slender 
shape  the  extrema  are  provided  by  the  choice  e$  = 0. 


Whereas  the  blunt  self-preserving  laminar  and  turbulent  shapes  described 
by  equations  (9a, b)  and  (15a,  a,  c)  essentially  coincide  with  those  reported 
In  References  1 and  2,  the  slender  turbulent  shape  of  equations  (16a, b,c)  repre- 
sents a novel  result  which  permits  rationalization  of  the  experimentally  ob- 
served sub-critical  turbulent  shapes.  Inspection  of  Figures  1,  2 and  4, 
with  attention  to  the  indicated  location  of  transition,  suggests  that  the 
shapes  of  interest  result  from  the  combination  of  a blunt  laminar  solution 
descriptive  of  the  nose  tip  region  together  with  a slender  turbulent  solu- 
tion descriptive  of  the  concave  contour  region.  The  assumption  of  sudden 
laminar-turbulent  transition  permits  direct  joining  of  the  two  solutions. 
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The  requirement  of  equal  recession  rate  at  the  joining  point,  together 
with  the  plausible  assumption  that  the  turbulent  boundary  layer  has 
negligible  thickness  (i.e.,e  = w/2)  at  that  point,  readily  provide  the  rela- 


tion between  the  length  scales  R^  and  of  the  two  solutions,  viz. 


Ct/Rw = {§ifr  [4(j  + ~ i)/n:lA>s  ( 


poU.  RN  ji*5 


*{or C4(j  + 1)i(Y  -i)/r]_1/Y.BVst)  Re0'" 

in  terms  of  the  conventionally  defined  transition  Reynolds  number 
Re.  = (p„  U s./y  ) . based  on  arc  length  s.  to  the  joining  point.  Slight 
generalizations  of  (17)  arise  with  different  estimates  of  the  turbulent 
boundary  layer  thickness,  or  equivalently  the  value  of  e for  the  slender  tur- 


rVS 


07) 


,1-5 


bulent  shape,  at  the  joining  point.  In  all  cases,  however,  the  geometry  of  the 


overall  configuration  becomes  controlled  by  the  magnitude  of  the  transition 
Reynolds  number  Re^,  which  sets  the  length  scale  ratio  C^/R^  [see  equation 
(17)].  An  example  of  the  transition  Reynolds  number  influence  upon  nose 
geometry  is  vividly  displayed  in  Figure  5;  the  family  of  combined  blunt 
laminar/slender  turbulent  solutions  is  characterized  by  the  foilwring  joining 
point  conditions:  e = et  on  the  laminar  side,  e * w/j  on  the  turbulent  side, 
parametrically  varying  magnitude  of  Re^.  It  may  be  noted  that  configurations 
of  Increasing  overall  slenderness  are  obtained  as  Ret  increases.  Yet  one 
length  scale  (either  Ct  or  RN),  and  the  recession  rate,  remain  underternined 
las  was  the  case  for  the  fully  laminar  and  fully  turbulent  solutions)  until 
the  value  of  e i 0 is  prescribed  at  the  station  where  the  overall  nose 
joins  with  the  afterbody. 


The  straightforward  requirement  that  the  nose  cap  the  given  afterbody  plays 
the  rather  unexpected  role  of  defining  the  Reynolds  number  range  over  which 
combined  blunt  laminar/slender  turbulent  solutions  can  exist.  As  equations 
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(9a)  and  (1 6b)  show,  the  ordinate  of  the  nose/afterbody  functions  is 

at  best 

(r,/RN)-  [r(et)/RN]  + (3/4)(Ct/RN) 


r*  • {[r(et)/s(et)]  + (3/4)  (ct  /RN)[RN/s(et)]}  (li 

(U«/v0) 

Once  Re^  Is  given,  the  right  hand  side  of  (18a)  and,  therefore,  the  nose 
diameter  are  uniquely  determined.  The  value  of  r$  so  defined  must  equal 
Of  exceed  the  radius  of  the  afterbody;  otherwise,  the  combined  laminar/ 
turbulent  solution  cannot  exist.  In  view  of  (T8a)  the  range  of  Reynolds 
numbers  Re$  * (Ujrs/vQ)  over  which  the  requirement  is  satisfied  depends 
on  Ret  according  to 

|r(V  1 Re«Re.  1 /r*V  + 3/4(Ct)[RN  ]\ Re  ,,, 

The  magnitudes  of  lower  (say  Re$l)  and  upper  (say  Res?  ) bounds  in  (18b)  depend 
on  j,  r and  Re^  fey  way  ©f  equations  (9a, b)  and  (17);  typical  values,  for 
4 » 1,  T »(7/5)  and  selected  values  of  Re^,  are  tabulated  below 


Ret 

!!ii 

*«s2 

rs/r(9,) 

l©4 

.921  x IQ4 

3.252  x 10“ 

3.531 

5 * ie4 

* i®4 

1.349  x 106 

29.30 

l©5 

.921  x 1®5 

7.462  x 106 

81.03 

5 * 10* 

4.50$  * l®5 

4.125  x 108 

895.8 

1©« 

.921  X 10* 

2.332  x 109 

2532.0 

The  calculated  values  ®ff  r%//r  feecome  extremely  large  for  Re^  5 x 104, 
suggesting  that  development  Otf  the  associated  blunt-laminar/slender-turbu- 

teat  CWifi  gyrations  may  fee  prevented  in  practice  by  reasons  of  nose  tip  struc 

Rural  imtegrity.  One  may  eamjeetur*  that  the  blunt-turbulent  solution  of 

Figure  J would  develop  under  those  conditions. 


j 
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I 

I 
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This  view  is  in  qualitative  agreement  with  the  experiments  of  Reference  3, 
where  development  of  slender  concave  shapes  was  observed  only  over  a 
range  of  unit  Reynolds  numbers  whose  lower  and  upper  bounds  differed  by 
a factor  less  than  two. 

To  summarize,  the  consideratioBs  of  this  section  rationalize  the  development 
of  slender  concave  shapes  on  ablating  nose  tips.  However,  they  do  not 
explain  the  different  flow  behaviors  observed  with  subcr.itical  turbulent 
and  critical  turbulent  configurations.  The  presence  of  unsteady  shock 
layers  in  the  latter  case  readily  suggests  that  the  key  to  the  problem 
cannot  be  found  in  a weak  viscous-inviscid  intersection  viewpoint,  as 


adopted  in  this  section;  instead,  strong  interactions  must  be  allowed  for. 
The  broader  approach  is  discussed  in  the  following  section. 


III.  UNSTEADY  STRONG  INTERACTIONS  ON  ABLATING  NOSE  TIPS 


Oscillatory  shock  wave  motions .such  as  are  exhibited  by  the  critical 
turbulent  shapes,  have  been  observed  by  several  investigators  working 
with  ablating  bodies  (Reference  3),  non-ablating  spiked  and  concave  conic 
bodies  (References,  S,  6 7,  8)  as  well  as  effectively  concave  conic  bod- 
ies as  are  encountered,  for  example,  with  massive  lateral  injection  on 
Slender  cones  (Reference  9).  The  typical  qualitative  sequence  of  events 
in  such  oscillatory  flows  is  illustrated  in  the  sketches  of  Figure  6,  which 
are  based  on  selected  shadowgraphs  for  the  M.  * 6.8  flow  around  a flat- 
headed cylinder  with  a nose  spike  having  length  equal  to  three  quarters  of 
the  body  diameters  (Ref.  5).  The  cycle  starts  with  the  boundary  layer  on 
the  spike  separating  near  the  nose  and  forming  an  approximately  conical 
dead  air  region,  which,  if  continued  to  the  face  of  the  body,  would  strike 
It  well  below  the  shoulder  (Figure  6a).  As  the  separated  layer  reattaches 
on  the  body  non-tangentially,  a shock  arises  at  the  reattachment  point  to 
turn  the  supersonic  flaw  just  outside  the  layer  up  the  nose  of  the  body. 
Depending  upon  the  negnitude  of  the  angle  6 (Figure  7)  that  the  separated 
layer  makes  with  the  tangent  to  the  body  at  the  point  of  Impingement,  the 
shock  Is  either  weak  or  strong.  In  the  latter  case  the  shock  stands  off 
from  the  point  of  impingement,  and  the  flow  behind  it  is  subsonic.  Inter- 
estingly enough,  the  onset  of  oscillatory  flow  is  observed  only  under 
these  conditions  (Reference  S). 

The  presence  of  a detached  shock  and  locally  subsonic  flow  at  reattachment 
has  a significant  implication  In  diagnosing  the  mechanism  of  the  oscillatory 
motion.  Dot  only  must  the  mass  flow  balance  be  preserved  in  the  dead  air 
region,  as  emphasized  by  previous  investigators  (Refs.  5,6,7),  but  also 
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the  choking  requirement  at  the  body  shoulder  must  be  satisfied  by  the 
external  inviscid  flow.  The  recognition  of  this  dual  requirement  dis- 
tinguishes the  present  analysis  from  those  of  the  previous  investigators. 

The  flow  patterns  sketched  in  Figure  7 clearly  indicate  that  the  geometry 
of,  and  the  pressure  ratio  across, the  detached  shock  are  controlled  by 
the  details  of  the  reattachment  process.  The  pressure  ratio  seeks  a magni- 
tude such  that  the  mass  of  air  reversed  at  reattachment  balances  the  mass  of 
air  scavenged  by  the  sheer  layer.  However,  the  process  ignores  the  choking 
requirement  at  the  shoulder.  Analysis  of  Figure  6a  shows  that  the  attached 
shock  lies  closer  to  the  shoulder  than  would  be  the  case  in  the  steady  flow 
over  a flat-headed  cylinder  without  the  spike  (Reference  10).  Under  these 
conditions  the  mass  flow  cannot  be  accomodated  through  the  sonic  section; 
a disturbance  must  then  propagate  forward  from  that  section  to  cause  a re- 
duction in  mass  flow  and  a simultaneous  outward  displacement  of  the  detached 
shock.  To  be  sure  the  forward  propagating  disturbance  increases  the  amount 
of  air  reversed  at  reattachment  with  consequent  bulging  of  the  dead  air 
region  and  the  associated  bow  shock  (Figures  6b,  c);  this  process,  however, 
constitutes  a concomitant  effect  rather  than  the  forcing  cause  of  the 
oscillation.  The  view  is  supported  by  the  following  discussion  of  experimental 
results  obtained  in  Reference  8 with  the  blunt  PANT-1  nose  configuration 
and  the  relatively  slender  proboTCiciid  (1-10-10)  configuration  (Figure  8)  at 
M»  - 3.0. 


The  flow  over  the  PANT-1  configuration  exhibits  a catostrophic  oscillation 
with  reduced  frequency  (Strouhal  number)  S = (fD/U.)  “ 0.227.  As  the  separ- 
ated free  shear  layer  has  a time-averaged  inclination  of  37°  relative  to  the 
axis  (Figure  9),  the  flow  deflection  at  reattachment  exceeds  the  maximum 
allowable  for  the  prevailing  Mach  number.  A strong  reattachment  shock 
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dtvelops,  which.  In  the  course  of  the  oscillation,  moves  between  most 
forward  and  most  rearward  loci  almost  parallel  to  each  other  (Figure  9).  The 
beck-and-forth  motion  covers  the  relatively  small  distance  of  0.2  0,  less  than 
the  estimated  detachment  distance  of  0.270  (measured  from  the  shoulder)  for 
•n  80*  half-angle  cone.  The  time  averaged  shock  location  (Identified  as 
Pitot  shock  In  Figure  9)  Is  consistent  with  that  estimated  for  steady  M.  = 3.0 
flow  over  the  80*  half  angle  cone.  The  time  averaged  pressure  distribution  over 
the  portion  of  the  nose  Included  between  the  average  reattachment  point  and  the 
shoulder  also  agrees  with  the  same. steady  flow  estimate.  Thus,  the  time-aver- 
aged results  are  Independent  of  the  detailed  mechanism  of  bubble  growth  and 
collapse.  Instead,  an  unsteady  choking  process  Is  Indicated  as  underlying 
the  oscillation. 

The  choking  Instability  viewpoint  is  further  supported  by  the  particular  magni- 
tude of  the  reduced  frequency  S.  The  associated  period  t * f-‘  (4.4D/U-)  = 
(D/0.23U.)  is  approximately  equal  to  twice  the  average  transit  time  of  a fluid 
element  across  the  inviseid  hypersonic  flow  over  the  nose;  thus, the  shock  oscilla- 
tion Is  geared  with  the  oscillations  in  site  of,  and  mass  flow  through,  the 
sonic  station.  In  line  with  the  transit  time  scaling  under  hypersonic  condi- 
tions (vis  N ^ 3 for  the  blunt  configurations  considered  here)  the  period  of 
oscillation,  when  present,  is  insensitive  to  details  of  the  nose  shape  and  to 
the  magnitude  of  the  free  stream  Kach  number  and  Reynolds  number  (References 
5,6,7  and  8).  However,  the  geometry  and  flow  parameters  influence  the  sta- 
bility boundaries  for  the  oscillation  (Refs.  5 and  7).  In  brief,  the  shape  and 
site  of  the  dead  air  region  are  controlled  by  the  flow  at  the  reattachmemt 
point.  The  very  low  velocities  in  the  region  endow  it  with  a fast  response 
tine.  Accordingly,  the  shear  layer  tends  to  move  so  that  standard  separation 
relations  (e.g.  Reference  6 for  laminar  flow  situations  and  Reference  11  for 
turbulent  flow  situations)  as  well  as  the  Chapman  reattachmant  condition  are 
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satisfied.  The  strong  reattachment  shock  generated  in  this  process  violates 
the  choking  requirement,  and  wave  disturbances  ensue  in  the  inviscid 
portion  of  the  flow.  The  waves  upset  the  equilibrium  in  the  dead  air  region 
and  cause  It  to  rapidly  respond  by  swelling  (shrinking)  as  flow  reversal  at 
the  reattachment  point  increases  (decreases).  The  sequence  of  events 
depicted  In  Figure  6 for  a spiked  body  then  follows.  As  air  flows  into  it,  the 
recirculation  region  begins  to  enlarge  causing  a nearly  normal  shock  at  the 
spike  tip  (Fig.  6b)  to  replace  the  initial  conical  shock  wave  (Fig.  6a). 

The  shock  wave  grows  from  the  spike  tip  as  more  air  is  fed  into  the  recircu- 
lation region  which  assumes  a blunt  shape  (Figure  6c).  The  associated  shock 
pattern  resembles  that  produced  by  a jet  of  air  from  the  body  injected  against 
the  main  stream  (Fig.  6f).  Simultaneously  , the  mass  flow  through  the  shoulder 
throat  decreases  and  the  reattachment  point  moves  up  the  face  of  the  body. 

As  the  choking  requirement  is  accommodated,  the  disturbance  from  the  sonic 
Station  changes  sign  causing  air  to  be  drained  from  the  recirculation  region. 
The  region  then  begins  to  collapse  and  the  shock  wave,  which  was  at  the  spike 
ttp»  moves  toward  the  body  (Fig.  6d).  In  the  process  separation  again  occurs 
4t  the  spike  tip  with  a conical  shock  formed  (Fig.  6e).  When  the  strong  shock 
returns  sufficiently  close  to  the  body  (as  in  Fig.  6a)  the  cycle  begins  again. 

The  question  as  to  whether  the  oscillation  is  due  to  unsteady  choking  or  an 
inherent  inviscid  instability  of  the  mean  flow  profiles  upstream  of  the  strong 
Shock  retains  to  be  examined.  High  frequency  oscillations  have  been  observed 
in  experiments  and  calculations  for  the  flow  field  around  the  proboscloid 
(t~lQrO(2)  configuration  (Fig.  8,  Refs.  8 and  12)  at  widely  different  Reynolds 
Mntoers(5  * IQ4  < ReQ  < 1.5  x 105  in  the  experiments  of  Reference  8 versus 
• A.8  x 103  in  the  calculations  of  Reference  121, different  free  stream 
Rich  *«mfe«rs  (Mw  > 3 in  the  experiments  and  M»  = 10  in  the  calculations) 
vut  t®  th*  presence  of  quite  different  flow  patterns  (separated  flow  in  the 
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experiments  versus  attached  flow  In  the  calculations).  These  oscillations, 
however,  are  quite  distinct  from  those  observed  with  the  PANT-1  configuration. 
The  reduced  frequencies  S » (fO/U«J  are  quite  large(4.5  £ S £ 9 for  the 
experiments  and  S * 8.2  for  the  calculated  fundamental  frequency);  shock 
wave  motions,  if  present  at  all,  are  quite  modest  (Fig.  10).  Measurable 
fluctuations  are  confined  to  the  portion  of  the  inviscid  flow  where  body  and 
shock  wave  converge  to  provide  a channelling  effect  and,  therefore,  spatial 
amplification.  Oscillations  with  the  noted  characteristics  cannot  be  asso- 
ciated with  a choking  instability  because:  a)  the  reattachment  shock  is  weak 

so  that  no  sonic  section  arises  downstream  of  it;  b)  the  periods  of  oscilla- 
tion are  smaller  than  the  average  transit  time  of  a fluid  element  across  the 
Inviscid  flow  region  where  measurable  oscillations  prevail  (Fig.  10).  In 
view  of  the  modest  sensitivity  of  the  reduced  frequency  to  Reynolds  rnanber, 
the  driving  mechanism  must  then  reside  in  an  inviscid  instability  connected 
with  the  presence  of  an  inflection  point  in  the  streamwise  velocity  profiles 
over  the  concave  portion  of  the  nose.  Such  an  inflection  clearly  exists  in 
the  separated  flows  encountered  in  the  experiments;  its  presence  in  the 
attached  flow  of  the  calculation  has  been  pointed  out  in  Reference  12.  Inter- 
estingly enough,  the  inflection  point  and  the  oscillation  disappear  if  calcu- 
lations are  carried  out  in  the  inviscid  limit  (Ref.  12). 

To  summarize,  oscillations  can  be  expected  in  flows  around  subcritical  turbulent 
as  well  as  supercritical  turbulent  ablating  nose  shapes.  The  former  are 
associated  with  relatively  slender  configuration  where  the  attachment  of  the 
separated  flow,  if  present,  is  accompanied  by  a weak  shock;  the  latter  are 
associated  with  blunt  configurations,  where  separation  always  arises  and 
reattachment  Is  accompanied  by  a strong  detached  shock.  In  all  cases,  viscous 
flow  processes  provide  the  qualitative  flow  pattern  features  which  lead  to 
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Instability,  viz.  Inflected  velocity  profiles  in  subcritical  situations  and 
strong  reattachment  shocks  in  critical  situations.  However,  driving  and  or 
amplifying  mechanisms  for  the  oscillation  are  of  an  inviscid  nature.  Thus, 
analysis  of  the  flows  in  question  requires  a strong  interaction  coupling 
of  viscous  and  inviscid  flow  regions;  however,  without  resort  to  extremely 
careful  description  of  the  transport  processes.  In  view  of  the  complex  flow 
fields  involved,  the  analysis  is  necessarily  numerical  in  nature.  An  out- 
line of  a relevant  approach  is  presented  in  the  following  section. 

In  connection  with  applications  to  flight  vehicles  it  must  be  noted  that  the 
solution  of  the  oscillatory  flow  field  problem  does  not  avoid  the  operational 
difficulties  associated  with  their  occurrence.  The  key  resides  in  nose  tip 
material  selection.  In  fact,  the  results  of  Section  II  indicate  that  develop- 
ment of  slender  concave  shapes  is  most  probable  in  flows  characterized  by  a 
low  transition  Reynolds  number,  i.e. , flows  over  nose-tips  whose  materials 
develop  critical  roughness  upon  ablation.  According  to  Figure  5 the  lower 
the  transition  Reynolds  number  the  blunter  the  shape  which  can  be  expected 
at  a given  flight  Reynolds  number  (based  on  afterbody  diameter)  and,  there- 
fore, the  more  likely  the  onset  of  critical  regime.  Conversely,  for  fixed 
transition  Reynolds  number  and  afterbody  size,  evolution  of  the  nose  toward 
a blunter  shape  and  a critical  regime  can  be  expected  as  the  unit  Reynolds 
number  increases.  A first  cut  determination  of  material's  adequacy  in  avoid- 
ing onset  of  critical  conditions,  for  given  trajectory  and  initial  nose  config- 
uration/size, can  probably  be  carried  out  by  means  of  simplified  heat  trans- 
fer analyses  (such  as  those  of  Section  II)  in  conjunction  with  available 
transient  ablation  codes. 
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IV-  COMMENTS  ON  NUMERICAL  FLOW  FIELD  MODELS 

The  discussion  of  Section  III  as  to  the  invfscid  nature  of  the  oscilla- 
tion mechanism  under  subcritical  and  critical  conditions  bears  directly 
on  the  selection  of  numerical  models  descriptive  of  the  considered  flows. 

As  pointed  out  in  the  discussion,  the  numerical  model  must  realistically 
describe  inflections  in  streamwise  velocity  profiles  in  the  subcritical 
case, and  formation  of  strong  unsteady  reattachment  shock  waves  in  the 
critical  case.  At  first  sight,  these  features  seem  to  require  use  of  a 
viscous  flow  model,  i.e. , the  numerical  integration  of  Navier  Stokes  equa- 
tions. However,  the  brute  force  approach  proves  to  be  impractical  even 
with  the  most  advanced  computer  facilities  because:  a)  the  flows  of  inter- 
est are  characterized  by  high  Reynolds  numbers  and,  therefore,  by  confinement 
of  important  viscous  transports  within  regions  having  thickness  much  smaller 
than  the  typical  nose  dimension;  b)  fully  three-dimensional  unsteady  flows 
must  be  considered,  as  distinct:  azimuthal  modulations  of  the  oscillation 
arise,  at  least  in  the  critical  case  (Ref.  8).  Thus,  the  Navier- Stokes 
attack  would  entail  a four  dimensional  problem  over  a domain  which  involves 
an  inordinately  large  number  of  mesh  points.  For  these  reasons  the  study 
of  Reference  12  was  restricted  to  axisymmetric  situations  and  low  Reynolds 
numbers;  the  associated  attached  flow  field  may  not  exhibit  the  interaction 
leading  to  the  critical  oscillation. 

The  viewpoint  of  this  report  is  that  extensions  of  available  inviscid  flow 
models  (e.g. , Refs.  13  and  14)  provide  the  key  to  a realistic  solution  of 
the  problem.  The  required  flow  field  features  due  to  viscous  processes  and 
interactions  are  embodied  in  the  calculation  by  Introducing  the  separating 
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boundary  layer  as  a source  of  vorticity,  which  is  convected  along  the 
boundary  of  the  recirculation  region  and  drained,  consistent  with 
mass  conservation  considerations,  by  a singularity  located  at  the  point 
of  reattachment.  Within  this  framework  transports  are  neglected  within 
the  recirculation  region.  The  overall  structure  of  the  flow  pattern 
is  attributed  to  the  inviscid  interaction  between  the  recirculating  flow 
and  the  external  flow,  which  includes  the  strongly  rotational  portion 
originating  from  the  separated  boundary  layer.  This  portion,  which 
retains  the  requited  features  of  inflectional  velocity  profiles  in  the  sub- 
critical  case  and  interactive  shock  formation  in  the  critical  case,  can 
, be  resolved  with  a reasonable  mesh  size,  particularly  in  the  important  region 

of  reattachment  where  the  lateral  scale  of  of  the  strongly  vortical  flow  is 
enhanced  by  the  large  local  adverse  pressure  gradients..  Empirical  justifi- 
cation for  the  view  resides  in  the  past,  reasonably  successful,  use  of  the 
model  for  studies  of  steady  near  wake  flows  past  slender  two-dimensional 
bodies. 
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V.  CONCLUDING  REMARKS 

The  development  and  aerodynamics  of  slender  concave  configurations 
on  ablating  nose  tips  have  been  analyzed.  Rationalization  of  these 
shapes  has  been  obtained  in  terms  of  novel  self-preserving  (equilibrium) 
solutions  of  the  nose  recession  problem  in  the  presence  of  transition 

from  laminar  to  turbulent  flow  regime  on  the  nose.  Based  upon  a study 

> , 

of  these  solutions  and  their  domains  of  existence,  the  development  of 
subcritical  and  critical  configurations  has  been  associated  with  the 
relative  magnitudes  of  the  transition  Reynolds  number  (based  on  wetted 
length  to  transition)  and  the  flight  Reynold  number  (based  on  body 
frustrum  diameter).  In  that  connection  an  important  role  of  surface 
naterial  and  its  roughness  has  been  noted. 

Upon  a ree*amimatien  of  available  experimental  data  the  distinction 
between  subcritical  and  critical  conditions  has  been  associated  with 
the  geometrical  characteristics  of  the  nose  sfiape;  specifically,  critical 
behaviour  t i.e.  large  shock  and  flow  field  oscillations)  arises  only  with 
blunt  configurations  characterized  by  boundary  layer  separation  near  the 
tip  followed;  by  reattachment  on  nose  portions  which  require  strong  re- 
attachment  shocks.  The  presence  of  distinct  oscillatory  modes  has  been 
noted  for  the  subcritical  and  critical  cases,  and  the  distinct  underlying 
inviseid  instability  mechanisms  have  been  pointed  out.  On  that  basis 
the  adequacy  of  numerical  analysis  by  extension  of  available  inviscid  flow 
field  models  has  be®»  suggested, and  the  nature  of  the  required  extensions 
has  been,  outlined. 
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Figure  6.  (a,b,c,d,e)  Schematic  diagrams  of  the  oscillatory  shock 

pattern  in  the  = 6.8  flow  over  a spiked  flat-headed 

cylinder,  (f)  Shock  pattern  around  a flat-header  cylinder 
with  front  air  injection  against  the  main  stream. 

After  Reference  5. 


cylindrical 

frustrum 


Geometry  of  nose  tips  investigated  in  Reference  8 
(a)  PANT-1,  critical,  (b)  Proboscioid  "1-10-10", 
subcritical. 


Figure  8 


Figure  9.  Map  of  tha  oscillating  PANT-1  flow  field  observed  ir  Reference 
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Rationalization  of  slender  concave  shapes  on  ablating  nose  tips  is  obtained  in 
terms  of  composite  laminar/turbulent  self-preserving  solutions  to  the  nose  receS' 
sion  problem.  The  occurrence  of  subcritical  and  critical  (violently  oscillatory) 
flow  conditions  is  discussed  and  associated  with  specific  geometrical  character- 
istics of  the  ablating  nose  shape.  The  presence  of  distinct  oscillatory  modes  ii 
noted  for  both  the  subcritical  and  critical  cases,  and  the  distinct  underlying 
inviscid  instability  mechanisms  are  pointed  out.  On  that  basis  numerical  flow 
analyses  by  extension  of  available  unsteady  inviscid  flow  models  is  suggested, 
and  the  nature  of  the  required  extensions  is  outlined. 


